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Abstract. A rigorous way to obtain sharp bounds for Stokes constants is 
introduced and illustrated on a concrete problem arising in applications. 
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1. Introduction 

Stokes constants (multipliers) relate to the change in asymptotic behavior of a 
solution of a differential equation as the direction toward an irregular singularity 
is changed (cf. If the constants are nonzero, then the change in behavior of 

the solution is nontrivial and this fact plays a very important role in a number of 
problems. 

Many interesting results arc known for linear ordinary differential equations, see 
m I13L I15L [T^ and references therein; papers [501 use hypcrasymptotic methods 
to express Stokes constants as convergent series. 

Stokes multipliers have been evaluated in closed form for a wide class of integrable 
systems sometimes using difficult and subtle arguments, see e.g. [HI 1121 1141 1181 IB?] 
and references therein. Integrability however is non-generic, and is thought to play 
a crucial role in any explicit evaluation of Stokes constants. 

A theory of wide applicability of extended Borel summation, Borel plane singu- 
larities and their relation witgh Stokes phenomena and was introduced by Ecalle, 
|1(J| . For generic nonlinear systems complete asymptotic expansions (transseries) 
of solutions and their Borel summability, singularities in Borel plane and formulas 
linking them to nonlinear Stokes transitions are rigorously obtained in The 
paper j2| finds the link between the structure of singularities of solutions in the 
Borel plane and summation to the least term, as well as with the behavior of the 
coefficients of the asymptotic series in generic nonlinear systems. 

In applications it often only matters whether the Stokes constants are nonzero, 
while their exact value is not relevant. The results in [H] and |7j provide a rather 
straightforward way of obtaining rigorous and accurate estimates of the Stokes 
multipliers for a large class of linear or nonlinear differential systems; in principle 
any prescribed precision can be obtained, as well as the information that a constant 
does not vanish if such is the case. 

One goal of the present paper is to present such a method and at the same time 
complete an argument in the proof by Tanveer and Xie )21| for the nonexistence of 
steady fingers with width less than ^ when small nonzero surface tension is taken 
into account. Their argument relies on a conjectured nonzero value of a Stokes 
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constant of the differential equation 
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(1.1) 2v"-t+ — = 

V 

Eq. (jLlf) appears as an "inner-equation" arising in the context of steady Hcle-Shaw 
cell fingers (see also |2). It convenient to illustrate our general technique through 
this particular equation. We show that two Stokes constants for (|1.1|) are given by 

^3/2213/14 

= ' V(l/7)r(3/7) 

^ 7r3/22i3/i4 
^'■'^ ^^ = "^^ V(l/7)r(3/7) 

with 

(1.4) l<fo<l + — 

In particular it follows from H1.4(l that indeed Si and 5*2 are nonzero"'^. 

It will become apparent that the method introduced here applies to generic sys- 
tems of equations whose irregular singularity has rank one (this is the most frequent 
type of irregular singularities in applications) . It relies on a detailed relation estab- 
lished in U] between the Stokes constants and the behavior for large index of the 
coefficients of the asymptotic series solutions, followed by inductive proof of bounds 
on the solution of the recurrence relation defining them. 

2. Stokes constants and exponential asymptotics 
Consider a system of differential system of the form 

(2.5) y'= (^-A-iB^y + fo(a-)+g(x,y) 

where A = diagA, A = (Ai,...,A„), B = diag/3, /3 = (/3i, /3„), fo = Oix^"^), 
g = 0{x~^; |yp; |a;~^y|) (as \x\ — > oo). with fo analytic at cx) on a half-line d and g 
analytic at (oo,0) under nonresonance assumptions |S] (a slightly weaker condition 
than the linear independence of the eigenvalues Xj over the rationals). The system 
(|2.5() has then a rank one irregular singularity at infinity. 

2.1. Power series solutions and exponentially small terms. The general type 
of formal solutions of differential systems in the presence of irregular singular points 
was studied in detail by Fabry [n| and Cope |S| (see also 0, 0, HI)- For 
(assumed nonresonant) the general formal solution for large x has the form 

(2.6) y = yo+ C'^c-'^-^-x'^'^Sk 

keN"\0 

where yo a-nd Sj. are power series (generically divergent) and C = (Ci, C„) are 
free parameters. 

The general formal solutions H2.6|l can in fact be calculated algorithmically, in a 
way that that will be briefly sketched. 

The power series solution yo is unique. To determine terms beyond all orders 
of this series, a formal calculation is carried out to find formal solutions which are 



^Numerical calculation gives b = 1.1722 ■ 
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small perturbations to yo! for instance one substitutes y = yo + ^ in (|2.5|l . where 
S ^ x^"^ for all q E M+ (as a; ^ oo on some direction). Since yo satisfies ()2.5|l we 
formally get 



where C is a vector of free parameters. 

Order by order perturbation expansion around 12. 7f) first produces power se- 
ries multiplying the exponential H2.7|l and then smaller and smaller exponentials, 
eventually leading to (|2.6|) . 

2.2. Transseries solutions. From the point of view of correspondence of formal 
solutions to actual solutions it was recognized that in general only asymptotic ex- 
pressions of the form (|2.6|1 {transseries, as introduced by Ecalle [TP|) can be lifted 
to actual functions. 

Transseries and their correspondence with functions constitute the subject of 
exponential asymptotics, a field which developed substantially in the eighties with 
the work of Berry (hyperasymptotics), Ecalle (the theory of analyzable fimctions) 
and Kruskal (tower representations and nice functions). 

Let dg = e'^R+ be a direction in the complex x plane. A transseries solution 
along dg is, in our context, a formal solution (|2.6|1 whose terms are well ordered 
with respect to the relation: ^ as a; ^ oo, x £ dg. In particular, a formal solution 
(|2.6() is a transseries along dg if constants Cj = for all j such that e"^^^ is not 
going to zero for x E dg. 

2.3. Exact solutions associated to transseries solutions. For any direction 
dg there exists a one-to-one correspondence between transseries solutions along dg 
and actual solutions that go to zero on this direction. The correspondence is nat- 
ural, constructive and compatible with all operations with functions (respectively, 
transseries) |10| . 

Actual solutions y corresponding to a transseries (|2.6|) on a direction dg have 
the same classical asymptotic expansion yg as x oo, x E dg. The constants Cj 
multiply small exponentials in the transasymptotic expansion of the solutions y, 
and are beyond all orders of the power series yo ■ They therefore cannot be defined 
using the classical Poincare definition of asymptoticity. 

For nonresonant systems with a rank 1 irregular singular point this correspon- 
dence was established in H] and [Hj. It was shown that the series yo is Borel 
summable (in a generalized sense); its inverse Laplace transform Yo(p) is analytic 
at p = and is (generically) singular at an array of points (which determine the 
Stokes directions). If we denote 





(5~e 



C 



(2.8) 




n=2 



then we have 



(2.9) 
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The other power series Sk (|k| > 0) m a transseries solution (|2.6(l are also Borel 
summable in the generalized sense. It is then shown in |H] that the series of functions 
obtained from a transseries solution by replacing the power series with their Borel 
sums is convergent for x large enough on appropriate directions. The solution 
thus obtained is asymptotic to the power series yo for large x on the direction on 
which summation was performed. The correspondence obtained between transseries 
solutions and actual solutions is one-to-one and is compatible with all operations. 

2.4. Stokes phenomenon. Consider a solution y that goes to zero for x oo on 
a direction dg. To associate a transseries y to y along dg means in fact to specify 
the parameters Cj, hence Cj may be different for different values of 6. 

It turns out. however, that the Cj, as functions of 6, are picccwise constant; the 
directions dg at which one of the Cj has a jump discontinuity are called Stokes 
directions (see [HI for more details). 

The way Stokes multipliers S are related to the classical asymptotic behavior of 
the solutions is given in the following Proposition of [Hj: 

Proposition 1. Consider eq. \2.5]) under the assumptions given. Assume Ai is an 
eigenvalue of least modulus. Without loss of generality we can assume that Ai = 1 
and ^(3j < for all j (these inequalities can be arranged in the normalization 
process, ^). Let 7* he two paths in the right half plane, near the positive/ negative 
imaginary axis such that \x^^'^^^e^^^^\ as x ^ 00 along 7*. Consider the 

solutions y of Y2.5}) which are small in any proper subsector of the right half plane. 
Then, along 7* we have, for some C, 

(2.10) y = (C ± i^Oeix-^i+^e-^^i + oieix-^'+^e-^^') 

for large x along 7^ (where ej is the j*'' unit vector Bj — (0, 0,__^^1^ , 0, ...,0)^. 

Stokes constants relate to the Maclaurin series of Yg in the following way ([H] 
and 13^): 

Theorem 2 (0). Under the same assumptions as in Proposition^ we have 

(2-11) Yr(0)= E ||i(;#^(^.e.+h,«) 

j;|A,l=i 

where Yq{p) is the generalized inverse Laplace transform of yo (see f2.8\) . 1^2. y\) ). 
4>j = arg'Aj (ordered increasingly, starting with Ai = \, ~ 0), hj(r) ^ r^^Yij-^ 
for large r. 

3. Main Results 

NormaUzation. To apply Theorem 13 to solutions of (|l.l|l . the equation has to 
be shown to be amenable to the normal form H2.5|l . The substitution 

(3.12) v{t) = t-^'^ ( 1 + u{x) ) where a; = y t^/'^ 

transforms (|1.1|) to 

„ 1 1 , 12 1 iu^ + 2u^ 12 1 

3.13 u ~u-\ u u -, T 

^ ' 7 X 49 x^ 2(1 -h u)2 49 x^ 



^There is a typo in formula (2.1) of |7|: /3' should read f3 (as obtained in its proof). 
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Substitution (|3.12(l is natural and a general procedure for finding normalizing sub- 
stitutions was described in 

Equation can be written as a system 

1 \ f \ 1 f 

7x / V 49 



where 

f ui \ , , , 12 1 3wf + 2ul 

U = , U) = ^ Ui r-!^ r^r 

V W2 y ' ^ ' ^ 49 X2 ' 2(1 

The dominant linear part of H3.14|l is diagonalized by substituting 



^ 1 + ^ 

142; ' 14x 



(3.15) u(x) = S{x)y{x) with S{x) = 

which gives the normal form H2.5|) with n = 2 and 

^'■''^ ^=(i -i) ' ^^(t 4 



14 

(3.17) fo(x) = ^fo ,withfo = ^(^ _\ 

(3.18) g(--y)-S^("/ -')y + iM-,^y)(-' 

Equation H2.5|l with (|3.1t)|) . H3.17|l . (|3.18|) has a rank 1 irregular singularity at a; = oo 
and it is written in normal form. The eigenvalues of the matrix A are A = (1, —1) 
and Pi=P2 = -1/14. 

The dominant power series in the transscries solution yo = (yo;i, j/0;2) is 

6 _2 , ,^,,87 _3 2028 _4 , ,,^57798 

yo-i+aix) = — X ^ + (-1)*^ X -{ X ^ + (-1) a; ^ + ... 

yu,i+^v J ^ ^ 343 ^ 2401 ^ ' 16807 

(ae{0,l}). 

Transseries solutions. For eq. H2.5|) . 13.16IIXTS)) the transseries solutions along 
dg with 9 e (— 7r/2,7r/2) must have C2 = (while Ci is arbitrary); transseries along 
directions with 6 E {tt/2,3tt/2) must have Ci = (and C2 is arbitrary). 
Stokes phenomena. The problem of physical interest depends on the solutions y 
which are classically asymptotic to the power series yo for |a;| 00 in the sector 
argx £ [0, 5tt/8] [2]. The fact that the asymptoticity is required on a large enough 
sector implies that this solution is unique and its transseries on any direction of 
argument in (0, Stt/S] are expansions H2.6|) with Ci = C2 = 0. 

Indeed, a transseries (|2.6ll for argx G [0,7r/2] must have C2 = 0. Since iR+ 
is not a Stokes direction however, C2 remains zero in all directions with arg.T G 
(7r/2, 57r/8]; but in these directions a transseries also has Ci = and since is 
not a Stokes direction, Ci — also for arga; g (0, 7r/2]. 

Since is a Stokes direction, Ci may become nonzero here, and its value is the 
Stokes constant 5*1. 

Relation for equation (ESI, (|3.16I3.18|I is 

(3.19) Y^"^) (0) = ^ ^'^^p ( 5iei + S2e-<'-+^^> + h(r) ' 
with h(r) ^ r~^ho for large r. 
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From H3.15|l we have u{x) = yi{x) + y2{x). Hence, using ()3.19|l . we get: 
(3.20) (0) = Y^}{Q)+Y^^{Q) ^ ^^'^^p [s, + S2C<'-+'A> + 0{r-^)) 

Though obvious, it is worth pointing out that the inverse Laplace transform Uip) 
of u{x) has the convergent series expansion at p = 

hence u''^^ {^) = if r is even and u''^^ (^i) = Ur+i > if r is odd. Then from (|OHl 
for r even it foUows that 

(3.22) S2 = -S'leTt*'^ 

Using for r = 2n - 1, and we get the formulas lO and ((0|) 

where 6 satisfies (|4.25ll . 

4. Estimating the constant h 

The constant h in (|4.24() below can be estimated within any prescribed accuracy 
by the procedure described in this section. For the problem at hand however, the 
inequalities (|1.4() obtained here are more than enough. 

Proposition 3. Any solution of equation J^S.H^) that goes to zero along R+ has the 
power series expansion 

(4.23) u{x) - ^ x-2"u2„ (x +00) 

ri>l 

where U2n > and 

(4.24) U2n = ^f(I77f(^"""^"r(2n + 1) (l + 0{n-^)) 
with b satisfying 

(4.25) 1 < & < 1.324 

Solutions of satisfying v{t) ~ t^^/^ as i ^ +00 have the asymptotic expan- 
sion 

00 

(4.26) -W-E^'^o = l 

where the c„ satisfy the recurrence relation: 
_ (7n-6)(77i-4) 



, , 1 (7n -7k- 6)(7n - 7k - 4) A 1 

(4.27) +2}^ 1 '-Cr^-k-l}_^C,Ck-^- -}_^CkCn 

" fe=l i=0 k=l 

Let dn be the sequence satisfying the recurrence 

(4.28) d„ = (Z!i^M!l^d,„_i (n>l), do = l 
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Clearly, d„ > and in fact 



Denote 



4 7 r(i)r(f) 



(4.30) 6„ = ^ 

From l|4.27() . 14.28|l the recurrence for 6„ is 

(4.31) bn = + Qn , bo = l 
where Qn = ^ Qn with 



n+ 1 (7n-7fc-6)(7n-7fc-4) ^ 

(4.32) = ^ 2 bn-k-lCn-k-lTk 

" fe=l 

(4.33) q-^_1_t; , Tfc = 2dfe6fe+T^ 



fe-i 

(4.34) U^Y. b^bk-^d^dk-^ 



Proposition 4. The sequence 6„ converges and its limit b satisfies the estimate 
Therefore 

(4.35) c„^r(2n+l)(^^) „-if_|^ (n ^ +«.) 

The proof of Proposition 0] relies on the following two Lemmas: 
Lemma 5. Let n>5. Assume there exist Ai,A2 > such that 
Ai <bk < A2 for all k with < k < n - 1 

Then 

(4.36) |Q„| < 4 
where 

(4.37) B = 0.6 + 0.0144 
The proof of Lemma |31 is given in M.ll 

Lemma 6. For Ai = 1 and A2 = 1.324 we have 

(4.38) Ai<bk< A2 for k^ 0,1, 2,.. .,7 
and 

B B 

(4.39) Ai + — <bk< A2- — where B = 1.0787, for all k > 8 

k k 
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The proof of Lemma El is given in M.2\ 
Proof of Proposition ^ This is an immediate consequence of Lemmas and 
Indeed, by Lcmma|H|we have 1 < &fc < 1.324 for all fc > 0. Then by Lcmma[Slwe 
have \Qn\ < Bn'^ for all n > 5, and by H4.31|l the sequence 6„ is Cauchy. The 
estimate (|4.25|) follows from Lemma El and the asymptotic behavior H4.35|l of c„ 
follows from H4.a0|l and H4.29|l . Relation (|4.35|l follows from (I4.29() and the Stirling 
formula. 

It only remains to prove Lemmas O and O 

4.1. Proof of Lemma [5l Note the following estimate for any > A^o > 3: 

1 didjv 1 1 

— > didN-i = 2 ; h — > didn-i 

dN ^ dN dN "H, 

^2^lj^ + (^_3)^Hj^< 1 ^(^)< 1 forA^>Aro>3 

dpf dM Jy Jy 
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49(Ar- |)(iV- I) 49^N - f)iN - f){N ~ f) 



where 

(4.40) E{N) ^ 
Therefore 

1 ^"^ 1 

(4.41) Dn^^^Y^ d,dN-^ < -TT^EiNo) ior N>No>3 

dN ^ N-^ 

It is easy to check that the estimate also holds for A^ = Aq = 2. To estimate 
(see (|4.34|l l note that since it was assumed that > Ai > for all fc < n — 1, then 
for Aq < fc < 71 — 1, we have, using H4.41|l . 

(4.42) <Tl< AlDkdk < AlE{No)^dk for Nq < k < n - 1 

k^ 

Let A^o = 5; we have E{b) < 0.24. For fc = 2, 3,4 to estimate further note 

that Dkdk = k'^Dk% < 0.22^. Then from we get 

(4.43) < Tfc < 0.24 Al^dk < 0.06 Ajdk for 2 < fc < n - 1 

fc^ 

It follows that (sec (gSSj)) 

(4.44) Q<Tk<adk, for 1 < fc < n - 1 , where a = 2A2 + 0.06 Aj 
For 1 < fc < n - 1 we have, using (|4.44() . H4.29|l 

(4.45) < bn-k~idn-k~iTk 

(7n-7fc-6)(77i-7fc-4) _^ 



(4.46) = 2A2dn-kTk < 2aA2dkdn-k 
and using (|4.41() (see l|333l) 

^, aA2 ''^ , , . aA2E(No) 0.24ayl2 

(4.47) 0<Q+<-^J2 dkdn-k = «^2i?n < 'J " < -—r^ 

k—1 
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for n > 5. The term (see ()4.3H|I ^ is estimated similarly, using H4.43|l : 

(«8) „<<3.;<4l|P<5:lidI,„,„>5 

Then (|4.47|) . H4.48|l imphes H4.36|l which proves Lemma [S] 
4.2. Proof of Lemma 1^ A direct calculation yields: 

169 743 426573 

bi = l , 62 = = 1-061 • • • , 63 = = 1.092 • • • , 64 = = 1.113 • • • 

'160 ' 680 ' 382976 

_ 71300607 _ -, -.26 b - 1406520669011 _ 

^ ~ 63289600 ~ ' '"' ^ ~ 1239463526400 ~ ' 

135335882622883 , , 6575066918153233021 

&7 = = 1 140 • • • ba = 1 144 • • • 

118668949344000 5744440195153920000 

Then (|4.38|1 holds, and also H4.39|l is true for n = 8. Estimate (|4.39() is shown by 
induction. 

Let n > 9. Assuming H4.39|l for all k with 8 < k < n — 1, then in particular 
Ai < bk < A2 for all k with < A: < n - 1 hence holds. Using l|OT|l in 

for fc = n — 1 we get 

13 13 

Al + +Qn<bn< A2 



n — 1 71 — 1 

which in view of (|4.36() implies ()4.39() for k ~ n. Lemma El is proved. 

Note 7. Substantially sharper estimates of b can he obtained using more terms in 
the expansion of'Yl^\o) that can be easily obtained from [Jj. 
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